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ABSTRACT: In this paper we formulate recommendations to decrease mismatch-induced stress in a 


heterostructure during microwave annealing of radiation defects, which were generated during 
radiation processing. Also we considered an analytical approach to analyze mass and heat transfer 
during the considered process. The approach gives a possibility to take into account changes in the 
parameters of processes in space and time, as well as the nonlinearity of the considered processes. 
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INTRODUCTION 


At the present time different heterostructures 
are widely used for the manufacturing of solid- 
state electronic devices [1-5]. The lattice 
constants of the heterostructure layers have 
different values. In this situation, mismatch 
induced stress’ existing in layers’ of 
heterostructures. Different methods could be 
used to decrease the stress. One of them is 
selection of materials with smaller difference 
between lattice constants. Another one is using 
buffer layer between main layers of 
heterostructure with intermediate value of their 
lattice constants. In the framework of this paper 
we consider an alternative approach to 
decrease the value of mismatch-induced stress. 
In this study, the approach influence of 
radiation processing of materials’ of 
heterostructure on value of mismatch induced 
stress is considered. 


In this paper, we consider a_ two-layer 
heterostructure. The heterostructure consist of 
an epitaxial layer and a substrate (see Fig. 1). 
Let us also consider the radiation processing of 
the heterostructure with future microwave 
annealing of radiation defects, which were 
generated during the _ processing. An 
accompanying aim of this paper is introduction 
of an analytical approach for analysis mass and 
heat transfer. The approach should give a 
possibility to take into account changes in the 
parameters of processes in space and time, as 
well as the nonlinearity of the considered 


processes. 
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Fig. 1 Heterostructure containing an epitaxial 
layer and a substrate. This figure also shows 
the qualitative distribution of radiation defects 
in a semiconductor heterostructure 


METHOD OF SOLUTION 


In this paper the distributions of 
concentrations of radiation defects in space and 
time has been calculated by solution of the 
following system of equations [6-8] 
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Here T (x,y,z,t) is the spatio-temporal 
distribution of temperature; x, y and z are the 
spatial coordinated; t is time; p =/,V; I (x,,z,t) 
and V_ (x,y,z,t) are the spatio-temporal 
distributions of concentrations of interstitials 
and vacancies with equilibrium distributions /* 
uw V’, respectively; functions D,(x,y,z,T), 
Dy(xy,Z,T), Dis(Xy,Z,T), Dvs(xy,z,T) describe 
dependences of volumetric and surficial (under 
influence of mismatch-induced stress) diffusion 


coefficients of point radiation defects defects on 
spatial coordinate and temperature; squares of 
concentration of point defects /(x,,z,t) and 
V2(x,y,z,t) gives a possibility to take into account 
generation of their simplest complexes with 
concentrations ®, (x,y,z,t) and ®y (x,y,z,t) (see, 
for example, Ref. [7]); function kiv(x,y,z,T) 
describes dependence of parameter of point 
radiation defects recombination on_ spatial 
coordinate and temperature; the considered 
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functions k;,(x,y,z, T) and ky,(x,y,z,T) describes 
dependences of the parameters of simplest 
complexes of point radiation defects generation 
on spatial coordinate and temperature; is the 
atomic volume of dopant; Vs is the symbol of 


L, 
surficial gradient; (C(x y,zt)dz is the 
0 


surficial concentration of dopant on interface 
between layers of heterostructure (in this 


00,(x,y,zt)_ 0 
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situation we assume, that Z-axis is 
perpendicular to interface between layers of 
heterostructure); k is the Biltzmann constant. 
Distributions of concentrations of radiation 
defects in space and time has been calculated by 
solution of the following system of equations [6- 
8] 
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Here Do(X,y,Z,T), Daev(X,y,2Z,T), Deis(Xy,Z,T) and Devs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of complexes of radiation defects; k;(x,y,z,T) and ky(x,y,z,T) are the parameters of 
decay of complexes of radiation defects. 

Chemical potential in Eq.(1) could be determine by the following relation [8] 


LW =E(zZ)Qo; [uy(Xy,2,0) +uj(Xy,2,0)]/2. 
1 Ou, ce Ou, 
2\ Ox, Ox; 


tensor; uj uj; are the components u;(x,y,Z,t), u(xy,z,t) and u,(x,,z,t) of the displacement vector 


is the deformation 


Here E(z) is the Young modulus, oj is the stress tensor; u, = 


u (x,y, Zt ); X;, X; are the coordinate x, y, z. The previous equation could be transform to the 
following form 


Ox, Ox, 2 Ox, Ox. 
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1-26 (z) Ox, 2 
Here o is Poisson coefficient; €0 = (@s-ae1)/azx is the equilibrium temperature, which coincide 
the mismatch parameter; ds, ds, are lattice (for our case) with room temperature. 
distances of the substrate and the epitaxial Components of displacement vector could be 
layer; K is the modulus of uniform compression; obtained by solution of the following equations 
B is the coefficient of thermal expansion; T, is [9] 
P ae u,(x, y, zt) _ do , (x,y, z,t) n 00 ,, (x, y,z,t) n 00 .. (x, y,z,t) 
of Ox Oy Oz 
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(z) 5 = : + = + (3) 
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Staley 28 B(z)K (z)[T(x, y,z,t)—T, ], p (z) is the density of materials of 
Xi, 


heterostructure, 5, Is the Kronecker symbol. With account the relation for oj last system of equation 
could be written as 
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Conditions for the system of Eq. (3a) could be written in the form 
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In relation (3a), both the nonstationarity of the temperature field and its inhomogeneity are taken into 


account. The space-time temperature distribution can be described using the corresponding Fourier 
law [10] 
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Here A is the thermal conductivity coefficient. 
Value of the coefficient depend on properties of 
the heterostructure materials and on value of 
temperature. The dependence of the thermal 
conductivity coefficient on temperature in the 
considered interval of temperature could be 
approximated as follows: 
MXY,Z,T)=Aass(X,y,Z)[1 +p Ta®/T?(x, y,Z,t)] (see, 
for example, [10]).  c(T)=cass[1-0 exp(- 
T(x,y,z,t)/Ta)| is the heat capacity of the 
material; Ta is the Debye temperature [10]. In 
the most interesting interval of values the 
considered temperature is comparable with the 
Debye temperature or exceeds it. In this 


6 T,(x, y,z,t) 


ass oC t 


situation one can consider: c(T)*Cass. p(X,y,Z,t) is 
the bulk power density, which released in the 
heterostructure during annealing. At the first 
stage let us_ calculate spatio-temporal 
distribution of temperature. Let us _ solve 
equation (4) by method of averaging of 
functional corrections [11-14]. Framework this 
method we_ replace the _ spatio-temporal 
distribution of temperature on the right-hand 
side of the equation (4) by its not yet known 
avearge value Our. The replacement gives a 
possibility to obtain the equation for calculation 
of the first-order approximation of the 
considered temperature in the following form 


= p(x, y,z,t). (5) 


Integration of left and right sides of Eq. (5) on time t leads to the following result 
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Not yet known average value 017 could be determine by standard relation [11-13] 
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The second-order approximation of temperature could be obtained by the following replacement T 
(x,y,Z,t) Or 2+T1(x,y,z,t) in the right side of Eq. (4). Appropriate equation could be written as 
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Integration of the above equation on time gives a possibility to obtain equation for calculation of the 
second-order approximation of temperature in the following form 
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Average value of the second-order approximation of temperature ar could be determined by the 
following standard relation [11-13] 
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Solution of Eq. (10) (ie. average value Gr 2) 
depends on value of parameter @. Next let us 
calculate solutions of equations of system (3a). 
These solutions describe components of 
displacement vector. To determine the first- 
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using method of averaging functional 
corrections, we substitute the not yet known 
average values a; into the right side of the 
equations under consideration. This 
substitution leads to the following result 


order approximations of these components by 
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Integration of left and right sides of the above equations on time t leads to the following result 
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Approximations with the second and higher orders of components of displacement vector could be 
obtained by standard replacement of required functions in right sides of Eqs. (8) on the following sum 
ait+ui(x,y,z,t) [11-13]. The replacement leads to the following result 
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We calculate spatio-temporal distributions of concentrations of radiation defects by standard method 
of averaging of function corrections [11-13]. To use the approach we transform Eqs. (1) and (2) with 
account initial distributions of the considered concentrations, i.e. 
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Further let us replace considered functions in right sides of Eqs. (1a) and (2a) on their not yet known 
average values Qip. The replacement gives a possibility to obtain equations for the first-order 
approximations of the considered concentrations 
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Integration of left and right sides of the above equations (1b) and (2b) gives a possibility to obtain 
relations for the first-order approximations of required concentrations in the required form 
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We determine average values of the first-order approximations of required concentrations by using 
standard relations, which is analogous to relation (7) (see, for example, [11-13]). Substitution of 
relations (1c) and (2c) into relation (7) gives a possibility to obtained required average values in the 
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The approximations with the second and higher 
orders of the concentration of radiation defects 
were calculated framework of the standard 
iterative procedure of the method of averaging 
functional corrections [11-13]. Framework the 
procedure to find the approximations with the 
n-th order of the concentrations of radiation 
defects, we replace in the right sides of 
equations (1a) and (2a) the’ required 


distributions of concentrations / (x,,z,¢), 
V(x,y,z,t), Bi(x,y,z,t) and ®y(x,y,z,t) by the sum of 
the average value of the  n-th-order 
approximation and the n-1-th approximation of 
these functions, i.e. for the Sum Onp+P n-1(XJ,Z,t). 
After such a replacement, we obtain the 
relations for the second approximations of the 
required concentrations 
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Integration of left and right sides of the above equations (1d) and (2d) gives a possibility to obtain 
relations for the second-order approximations of required concentrations in the required form 
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Average values of second-order approximations relations (1e) and (3e) into relation (9) gives a 
of the considered functions were calculated by possibility to obtain following relations for the 
using relation, which is analogous to relation (9) required average values 0 2p 


(see, for example, [11-13]). Substitution of 


7 Oa,F +@°L,L,Lb, 
O2@;=0, A2z@v=0, Oey = (b+ EY —-4/Fr+ a, no ty *z _ Db, +E . 
Ab; b, Ab, 


Cy Oy Swen —Qoy (2 Swo1 ale Sry10 +OL,L,L )- Svv02 —Sivi 


x “yz 


Q,) = 5 5 
voi + oy pv00 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.4. 2431-2450 (2022) 2448 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


! 2 1 2 


Where b, = ar VV00 SELL Tw 00" 
xyz xy =z 
b, = (285 tp er i @L,L,L,) x 
a eS Oe On. 
OL LL OLL L800 2510 * Sat OL LL + rpg Si 
= Si) ne 
OLLL 
Sanat > SiS 
b, = @LLL. + Sit + C,)- (OL,L,L. ~ 28 wii + Svio) an OLLL. x 
OLLL +S,,.+2S8,,, +258 
2 (GL LT £98..¢5: 6 06 POLIT AS.) x 
OL,L,L. 
x a —S? 4 Cy 7 Sate 7 De r CS ci mp Die” aaah wai 
OL.L,L. @LEL OL.L,L. 
b= Car + Syyo2 + C,) x 
28 yo, +S, + OL LL, OLLL,+2S,,,+5 
x Sion = ait L — _ Sg = OL L - _ (284 + Si +OL,L,L,)- 
x yk x yz 
C, -S,,. —-S a 
x br). Se : ore L a (35; + 28 no + OL,L,L,)+ 2C Si Wai ~~ ane 
SS a) CH5 328 
b, = Ona ( ol L ne Dai : an L = (OL,L,L. PLS ig a Co Ie. el _ 
xy x 2 ¥ 2 
C,-S,..—S oa, oO aS 
vol . OLLL. a (0 LLL, + LD vais ae Sa) C, = OLLI" a OL.LL, 
SS oq £9 
= nao Cy = 6,017 S iv00 70,570 a Sywor — Sivir 


2 
b= by+0° 4-40 


a, a, 
Oa 
Pov tlwr +s" —r-\vr+s +r, 
a, 


3 2 
po lag ere 1 sp os 
24h: ee 8b: 


4 


Pankratov 


International Journal of Advanced Science and Engineering www.mahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.8 No.4. 2431-2450 (2022) 2449 


poe 
8b 


4 


2 32,3 
x| 405, -97 2 |-2 pr 
b,) 54b3 


4 


In this work, the required distributions of 
temperature and concentration of radiation 
defects, as well as the components of the 
displacement vector were calculated by using 
the second-order approximations framework 
the method of averaging functional corrections. 
The approximation is usually enough good 
approximation to make qualitative analysis and 
to obtain some quantitative results. All obtained 
results have been checked by comparison with 
results of numerical simulations. 


DISCUSSION 


In this section we analyzed relaxation of 
mismatch-induced stresses and spatio-temporal 
distributions of radiation defects. It was 
obtained, that it is practicably to choose energy 
of the implanted ions so, that maximum 
distribution of radiation defects (see Fig. 1) 
would be inside of the heterostructure after the 
interface between its layers. Taking into 
account that the value of the diffusion 
coefficient of interstitial atoms is usually larger, 
than value of diffusion coefficient of vacancies, 
with time, in the neighborhood of interface 
between the heterostructure layers, quantity of 
radiation vacancies will larger, than quantity of 
interstitial atoms. In this case, at €o>0, the value 
of mismatch-induced stresses decreases with a 
simultaneous decreasing of quantity of 
radiation vacancies. The dependences of the 
values of the displacement vector and the 
concentration of vacancies on the coordinate, 
which is perpendicular to the interface between 
the heterostructure layers, are shown on Figs. 2 
and 3 respectively. Probably, the decreasing of 
quantity of vacancies became more compact in 
the epitaxial layer under influence of mismatch- 
induced stresses. At the same time with 
increasing of density of the epitaxial layer 
decreasing of value of mismatch-induced 
stresses could be obtained, probably, by 
increasing of the density of atoms of the 
epitaxial layer. Analysis of distribution of 
temperature shows that it is advisable to choose 
the frequency of electromagnetic radiation so, 
that the thickness of the skin layer would be 
comparable to the width of the irradiated area 
by radiation processing. Decreasing of the 
considered frequency is not attracted an 
interest, since the energy of the electromagnetic 
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field will be spent on the area of the 
heterostructure, which was not damaged during 
radiation processing. Increasing of frequency 
leads to the fact that the deepest area damaged 
during radiation treatment will not be annealed. 


Fig. 2. Normalized dependences of the 
component of the displacement vector u, on 
the z coordinate for epitaxial layers before 
radiation processing (curve 1) and after 
radiation processing (curve 2) 
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Fig. 3. Normalized dependences of vacancy 

concentrations on the z _ coordinate in 

unstressed (curve 1) and stressed (curve 2) 

unprocessed materials 


CONCLUSION 


In this paper, we consider the effect of 
microwave annealing of radiation defects 
generated during radiation processing of a 
heterostructure with account relaxation of 
mismatch-induced stresses. Based on_ the 
analysis we formulate recommendations to 
decrease the considered stresses. An analytical 
approach for analyzing of mass and _ heat 
transfer is introduced. The approach gives a 
possibility to take into account changes of 
parameters of processes in space and time, as 
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well as the nonlinearity of the considered 
processes. 
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